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ABSTRACT 

An experiment 1S under way which will attempt to 
measure the neutron-neutron quasi-free scattering cross 
section in neutron-deuteron break-up to an accuracy of a few 
percent. This work is a theoretical study of the reaction as 
a means of determining the neutron-neutron effective range. 
Both impulse approximation and three-body models have been 
used. Some recommendations regarding the experimental 
geometry are made. The ingredients of a proper theoretical 
error analysis are enone 

The effective range can be deduced from the 
experimental cross section with a percentage uncertainty 
equal to that of the cross section. A further theoretical 


error of.-2% must be included. 
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Lee NTRODUGTION 

The nonexistence of a free neutron target has always 
made experimental investigation of the pure neutron-neutron 
(n=n) interaction difficult. Details of the low energy n-n 
interaction remain much less clear than those of the 
HMeLUELON PrOoLronnand Puoeon=protone@interactions.Nirhough 
important in their own right such details could supply 
notable tests of charge independence and charge symmetry. 

One obvious place to attempt to observe the n-n 
interaction is in the three-nucleon system, either in the 
bo@nd Statea@(triton) of @ninetibronsdeuterone (n=D) 
Scattering. The difficulty here, as in any interacting 
many-body system is that the off shell nature of the 
constituent two-body interactions hampers the observation of 
on-shell two-body scattering behaviour. 

Chew (@€h46550052)eseemsstoahaveubeenythestirste tovlook 
for two-nucleon information in the three-nucleon system, 
trying to deduce the free n-n and n-p cross sections from 
n-D inelastic scattering results using the impulse 
approximation. Since then this method has been widely used, 
and with reasonable success, to describe n-D elastic 
Scattering (Ko65 and references therein). Success with n-D 
and p-D inelastic scattering has been less (references in 
Du73) especially at low energies. It was not until practical 
formulations of the three-body problem (Fa60,61,63,65, 
We63,64, Lo64) and subsequent calculationsf that accurate 


These had to await a further breakthrough in high-speed 
computing. 
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determination of n-n scattering parameters became possible. 

n-n final-state interaction (FSI) experiments have been 
POUnGHSenSi)t1Ve Moe thesm=n Scatceringmlengtn (aye), whech 
has been determined to an accuracy of about 3% (Ku76, Wi79, 
Ss81). Using phase-equivalent potentialst Haftel (Ha73) and 
Harper (Hr72) have shown that the triton binding energy (E,) 
1S sensitive to the off-shell behaviour of nucleon-nucleon 
(N“N)Sainteractions. The Phillips line (Ph68) demonstrates a 
direct relationship between E, and the n-D doublet elastic 
peer tee length. Brayshaw (Br74) has further asserted that 
there is essentially no off-shell dependence of the 
three-nucleon system beyond that specified by E,. This 
rather strong statement has seen considerable debate 
RE Sy iho. the dhsay  Snetie IAs) 

Quasi-free scattering (QFS), in which one of the target 
particles remains at rest, iS particularly appealing for 
SunayingeLhewm—-Meimteraction. Caloulations and Iantuation 
both suggest that this two-body-like process is almost 
entirely insensitive to the off shell behaviour of the N-N 
interactions. Moreover a strong sensitivity of the n-n OFS 
cross section to the n-n effective range (r,,) has made this 
a popular place to measure r,z,. Determinations of r,, (table 
wp eocestypicalily much less accurate ethan Chose Otea.,, tne 
Mimering factors being experimental. (Fou this reason, sand 
because of the generally low sensitivity of n-n OFS to other 


tpotentials that agree on-shell but have different. off-shell 
behaviours 
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aspectsy; of the n=n interaction, analysis of theoretical 
errors has not been extensive. Comparing a measured n-n 
effective range with the Coulomb-corrected p-p value would 
provide a simple test of charge symmetry. The corrected n-p 
(singlet) value provides a similar test of charge 
independence. Cursory inspection of table I.1, however, 
reveals that considerable improvement in accuracy is 
required before any meaningful comparisons can be made. 
Experiments other than n-n OFS are in evidence in table 
I.1 but none appears to offer any advantages in accuracy. 
The results from FSI measurements are considered unreliable 
due to an inadequate model and the low sensitivity to r,, in 
this region (Gu80). Of more interest is the report of von 
Viueschmetwal= s(WiSO)ewho claim that the ratio of n-p to n—n 
OFS cross sections can be meaSured more accurately and Phy 
deduced from it with less theoretical uncertainty. While the 
first claim may be true, their experimental uncertainty of 
8% leaves much room for betterment by conventional methods 
(Ss78). The second claim is discussed in chapter IV. 
Mheslatestesurvey by Slaus (Ss/6,.80) points ocutea 
systemmatic discrepancy between effective range parameters 
deduced from n-n QFS and (z~,y) experiments (on D). Slaus et 
aleecSsS?) = propose that a three-nucleon) force, implicit in 
themnature of the Strong, interaction, 1s) present. in the 


former but not the latter where only two nucleons 
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participate. One might, therefore, regard the differénce 
between 7,, measured from two-and three-nucleon experiments 
as a measure of the presence a three-nucleon force. The 
dominant component of such a force would be a two-pion 
exchange tterme(Mk76) } "Slobodrian “(Sb82)) maintains thar; 
short of colliding neutron beams, the only reliable way to 
Study Wvhecn=n interaction Wwsiiviaicempar ison iofe"minnor 
processes" in which neutrons and protons are interchanged. 
He cites Gross et al. (Gr70) who obtained agreement with 


(m-,y7) results using *He(*He,‘*He)2p and *H(?H,‘*He)2n. The 
comparison method has been criticized (Ss71). 

Sauer and Walliser (Sr77) have found that the process 
of making electromagnetic corrections to p-p data is model 
dependent. They propose using the assumption of charge 
symmetry (i.e. equality of n-n and corrected p-p parameters) 
as a constraint on the off-shell behaviour of the nuclear 
PeOLreneialvelhe correction of n-p data does not appear to 
suffer from model dependence. 

Here at the Nuclear Research Centre an experiment is 
currently in progressft which will attempt to measure the n-n 
OFS cross section to within a few percent. This should 
Eiable aN determinacuonsot MawetOla similar previously 
umatcained) degree of accuracy. To maintain this level of 


accuracy, a more comprehensive and detailed theoretical 


understanding of the reaction is required. The purpose of 


+Details of the experimental setup can be found in section 
Mids GeanGeinuNucl camenesecanene Gentine sProgresce Reports 
1979-82. 
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this work is to provide that understanding. Of course the 
relationship between r,, and the OFS cross section must be 
well established. The sensitivity of the outcome to 
experimental geometry and resolution must be explored. 
Finally the various theoretical uncertainties that have 
previously been swamped by experimental error must be either 
eliminated or quantified. The Sasi 6: the investigation 15s 
dictated by the anticipated experimental uncertainty. 

Such detailed analysis is only possible with the latest 
three-body computer calculations. We have used a code 
supplied by Pal Doleschall (Do73) which we believe is the 
best available for the purpose. The results of its use are 
presented in chapter IV. Preceeding these (in chapter III) 
are the results of our own "poor man's" three-body 
calculations using the impulse approximation. Chapter II is 


a review of basic three-body theory, intended mainly for 


perspective. 
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Il. THREE-BODY THEORY 


a. Notation 

In three-body equations, the interacting particles are 
labelled 1,2 and 3. Quantities defined on a pairwise basis 
elem anetved by the index of the non-participating particle. 
Thus V,; indicates the potential between particles j and k 
Gime Ke cyClic@im This <dilsouss von), 

To simplify later expressions, Lovelace variables 
(Lo64) are used with all masses set equal to 1. Pairwise 


relative momentum is denoted by the vector quantity D} 
pe=Gis2) Gkr=k®) Dips? | 

where kK; is the laboratory momentum of particle i. en is 
proportional to the momentum of particle k relative to the 
(ty) eCentre-of-mass: 

Gua 1/2V3)( kK, +k; -2k,). Liga? 
A free state is labelled by its momenta, and in the 
three-body case can be expressed in four different but 
equivalent ways: 


lgi>|pir>(=|qipirJ=|kKikzks> I1.3 
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state 1S Simply expressed: 
E=p; 2+q; 2 bii4 


\ereDreSentsha® boundestaterormiimkieuniawe-body space. 


Accordingly, 
1s the same bound state in three-body space. 
Thesesymbolljagrandicatesuthe f=matrixafor®’ (j;k) 
scattering, while the arguments supplied indicate whether it 
lswanmatrixeelement) or! operatorgacting in two- or three-body 
Space. In two-body space 
TROD D7. fie =< Darien Ed bn > Lic 
where E is the two-body energy. In three-body space 
WCPO Pe Gi EJ=<G7 | <Dalt eon | qa 11.7 


where E is now the total three-body energy. Thus 
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Singularity structure, depends on whether it is used in two- 
or three-body space. The incorporation of other two-body 
operators into three-body space follows analogous lines. 

The complete potential acting in the three-body system 


is assumed to be the sum of pairwise potentials 

Var e Dye oS 
although there 1S no reason why the summation could not 
include a "three-body potential", involving all three 
particles simultaneously (Fu57, Ph66, chapter I). 


In analogy to the two-body case there is a Greenf 


function for the entire process 
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describe the pairwise interactions. The bound states |I> 
(equation II.5) are the eigenstates of H;. 

In analogy to the two-body t-matrix, there is assumed 
Eoebe wast - Metrix for 020 "scattering (0 denoting the state of 


three free particles): 


T=V+VGV | II.14 
=VtVGo iE 


where Go describes the propagation of the O state. 


b. Limitations of Two-Body Theory 


The Lippmann-Schwinger equation, 11.14 or 


Y=ot+GoVW IG hy ee 


belongs to the general class of Fredholm equations. Viewed 
as an operator equation, certain Statements about its 
solution can be madet provided the kernel operator GoV is 
Compact ime nUMer ical SOlutd ONS, sli! Ue smuUStd lyeconvc luted 
to a set of matrix equations for which compactness 
guarantees a finite number of equations. In practical terms, 
therefore, compactness is elevated to a requirement for 


SO lMIGLOMN. 


+Fredholm theorem. See for example Kreyszig (Kr78). 
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Figure 1). biPwo-=body scattering in two-body Space is 
a connected process. 


In the Hilbert space of quantum mechanics, compactness 
1s roughly equivalent to square integrability (boundedness 
Ment hemes chinitatenorm je eMore intirbtive by. attite tkernél tom 11.15 
1S compact for any connected process in which no particle or 
group of particles can be singled out as separate 
(noninteracting) from the remaining particles. 

Sinplesguwo-hbody Scattering, depicted in figure 11.1, 1s 
clearly a connected process, for which equation II.15 is 
soluble, as long as V satisfies certain basic 
requirements.f Strictly speaking, the kernel GoV for most 
potentials is compact everywhere but at bound state energy 
poles and on the right-hand energy cut. Nevertheless the 
problem can be dealt with rigorously by using a few 


*See Taylor (Ta72) for example. 
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Bigure li.2 Two-bedy Scattering in the presence of a 
Spectator particle is a disconnected process in 
three-body space. 


mathematical tricks. These are reviewed by Lovelace (Lo64). 
The situation in three-body space is more forbidding. 
Even the simplest process — two-body scattering in the 
presence of a Spectator particle — is obviously disconnected 
(figure 11.2). Mathematically, writing the two-body Green 
function in three-body space (following along the lines of 


equation II1.8) 
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introduces a $-function in the third coordinate which leads, 
upon integration, to a very awkward infinity. 
In more general fully-interacting three-body problems 


the above divergence arises again, and not just as an 


133 


isolated factor. Decomposing I1.14 using I1.9 and iterating, 
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reveals that there are disconnected terms (for i=j=...) at 
each order of the expansion (figure 11.3). The insolubility 
of one such problem is.discussed in more detail by Afnan and 
Thomas (Af77). One is forced to conclude that the 
established mathematical formulation of two-body scattering 


is unable to describe three-body processes. 


c. Three-body Equations 
Faddeev (Fa60,61,63,65) approached the mathematical 
problem of disconnectedness by proposing a different 


decomposition of the 070 t-matrix: 
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The 7' are found to satisfy "Faddeev-type” coupled 


equations: 
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Figure II.3 Expansion of the 070 t-matrix in simple 
two-body interactions leads to an infinite number of 
disconnected terms. (equation I1.17) 

tf Sie eee er eee a ei ea” 


eiqaiz ni a7 +280 a GeO eusiiag nolan J 
25 qadpend a@sart int i. a a al fous Ped 
: ei Se ,aerae me 

° - 


Ages ore = 
—_ 


Re) 


ot 


b=TtlaGelieri.? Tet 
=[ed 1 Gols Hla Gol 

+T iGol ;Gol ;+T ;GoT ;GoT 
eGo Gel ;AT -GolaGols 


+ 


with connected kernels. Associated with these are Green 


Funct rons 
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which satisfy similar equations: 


G'=(G,=-Go)tGol; (G'tG"). M2 


II1.20 suggests that G' (or T') describes the process in 
which tar toes ny and, kwhinteract diast ~i2hissiprocess gis 

clearly subdivided by equation 11.19 into connected and 
disconnected parts.ft For convenience the completely 


connected operators 


if tothe ails TENE, aes 


are then defined, satisfying 


ties 


PRecaumethe -o-funetion of equation 1158 which causes 7) °to 
De Esingula te. 


Ure any Te"? — 


>» .) - Sean; wo Ve, 42 


- 

. 
. _- + — 

7 a CT vr? 

~ - ~~ 
j i | YY a ° 
’ 
j - v a) % =| , 4 
2 

om. TI ;3 


TT 
2am 6 : 


bet l 4 
' i 
Re 289 7 2e@gure 
Ss i! 7 ‘ aA ~~ se if 24 U3 3See 
Ot Ly » a oe ' a | i ; } ‘Sy aT = iiwavue 


i > 
: 


beso 
ee 


Ke ee pt 


16 
Vale Go len lg tye) 


Intuitively, Faddeev has regrouped the multiple 
SCacrtering series for the O0+08 interaction such that the 
Printy Of Ralsconnegteaumterms. In equations telly areal 1 
Goimected into the first term (the two-body f-matricies) of 
Sguareon ee Toei icgurerm) 1-4.) 


For most calculations the 7' are further decomposed: 
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From the elements W;; which describe the connected parts of 
the 070 process, can be derived (Fa61,0s71) the amplitudes 
for the other possible processes, although the connections 
are not simply made. 

To make the derivation of amplitudes for the general 
1(i,j)Jt process more transparent, Lovelace (L064) chose to 
begin rather than end with these amplitudes. Depending on 
the use of incoming and outgoing scattering wave functions 


in the derivation, two operators are found: 


trefer to equation II.5 
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Figure I1.4 Faddeev has regrouped the multiple 
scattering series so that disconnected terms are 
more easily dealt with. (equation 11.19) 
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Lovelace showed that these satisfy Faddeev-type equations 


(the "Lovelace equations"): 
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which are compact and soluble but do require both potentials 
and sfalainicwes Bors ibhe) two; bod yjeunt enaetions-: White 
differing off-shell, both operators yield the same t-matrix 


for the on-shell process 
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as required. 
The problem of off-shell ambiguity was addressed by 
Alt, Grassberger and Sandhas (A167) who defined a single 


type of operator 


ion Gurls GS paul ona) 11.29 


which alsov describes the: on-shell three-body, process. In 


addition to being more symmetric, these operators satisfy 


Faddeev-type equations (the "AGS equations") which do not 
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involve explicit two-body potentials 
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The three-body breakup transition operator can be defined in 


EGrms of the bound state transition operators: 
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The impulse approximation used in chapter III can be 
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d. Separable Potential 

While three-body equations have been solved for a local 
potential) (Os67,K171)) the mumerical problem as, in general, 
very difficult.+ Even after angular momentum decomposition, 
which reduces the dimensions from six to two, the remaining 
integral equations are coupled and infinite in number. 

Most calculations, therefore, introduce separable 
partial wave two-body potentials (the partial wave index is 


suppressed): 
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Although this may seem rather restrictive, it is worth 
noting that any local potential can be approximated with 
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*There are mathematically sound theories (e.g. We64) but for 
Pract calecalculationsssimplit ying eassumptionc are required 
(see Sc74). 
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particular, the two-term R-type separable potential 
described in chapter IV has a level of agreement with 
experimental N-N phase shifts previously available only with 
non-separable potentials such as the Reid soft-core (Re68). 
In the presence of a separable potential, the See” 


partial wave t-matrix itself becomes separable: 
MC Deere y= 20) Cooter JOzep.e); iimsG 
PNeOpebhacor Nnotacron, 
T(E)=|g>r(E)<g| It-32 
where the conventions: Of Section A»still apply, i.e. 

Ge NE MWelep 5) ten. ten ey, Game corns) ep, PLE SSNs} 
MMmeeGee=bOOve Shace wt tOmSimMplity; moOtativon cles un ura 
variable will be omitted from equations in three-body space 
on the understanding that a g; should be present for each 
subscript i, with three-dimensional integration over 
repeated subscripts. 

In this notation, the "X-amplitude" for rearrangement 


Scare ning Tis ide fined: 


Xp=<Gn| Goal b Ui; Ge GE) guar ie Be 


ne tia 
sae Saute wa 


oval ia Seri V1 


Aste yine eat iava ‘aSetaties a 


; ~~ ' 4 tee 
tc ik i ae) we Lore { ep>! an 3 2G Baw © 3 o ~ : 
e = 4 
i> 5171427. * 
t * a 4 x Sa as se |= 
Gualov: ae ‘ te oe I 2-Ghd ° 
x 
4 4 oe so 
jiknstacad Bite ed Lf 
te oiik 
{ rmit« aj i - 4,4 ' 
7 
4 16707) F 
on 
>. j j } a = ; 
«4 {,t.~ @a 
’ 
a ' ¥ ; - { a] 
* 
f : 7 
Sted 4 E Aaa iit ened? 
7 ry é 
Ssunce ‘ al = ~ iI , Ss ; * on re * » 
7 
Ish5 201 | Bn Bs | au 6 HENS 7.4 J fn) ieee’ 2 
7eV0 WicaA i i GOMZe eels! ) a” ae | nTtoe 


acnlvoradpe 


sqaneenet bass. 20! ““ehu ) 2igqne-\"- siz. , chidetean ei oe 
5 7 j - " 7 


shanti ioB « 


» 


: 7) = 
; ‘ Ae ¥ 7 —— a 
- Fie ea i aj BS 7 
‘ 7 
aan 


Bs 


nse ren ngithe separable ife=mat rix fof W187) into 11.30 and 
taking matrix elements <g;|Go(E)...Go(E)|gi> yields a set of 


coupled mntegralcequatronsimin sone 'contunuoustvariable for 


Com: 
XC RGE)= 5 OV CAtElse eapemace Jan, Glade? 16 -a( oan 11. 40 
Vi; 18 an energy-dependent "effective potential": 
View G=407, Goll) loa. Tie] 
By taking matrix elements <po|...Go(EJ)|g;> of equation 


I1.31, analogous equations for the breakup case are obtained 


Xo j=Vo0 jl SOVTHZ Von TeXn j 11.42 

where 
Ko (B= ales! AGE) lg i> U1,<43 
View (ED =D lg oe Tei4.4 


After angular momentum decomposition (Do73,A167), 
equations II.40 become one-dimensional integral equations 
which can be solved by a variety of methods (some of which 


are described in section F). The breakup amplitudes Xo, can 


+See equation I1.49. 
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Al Choughitt hes tirmatrix ofmequat honsiliise6 and el .37 
could be expressed in any basis, the ideal (We63,64,Sc74) 
basis of bound state form factors is used to simplify 
calculations. In this representation, each form factor for a 
pairwise interaction can be related on shell to a bound 


State, Ole the parr: 


COa Gave GaP =< 0 | xq: Il .45 


louSea, DOLentlal will typically have tone form factor icr 
each bound state. 

The relation 11.45 can be used to attach some physical 
meaning to the mysterious X-amplitudes. From définition 
Peo 1G 1s Clear cthation ighnelik, fe) fis) ithe amplitude for 
scattering from a bound ‘state of (i,k) to a bound state of 


(Gia 

Xray eX Ui kee T1.46 
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Of course),,1if%a particlel pair hadhkmoresthan? oneabound 
state, X would have to be labelled accordingly. For 
Simplicity and because the deuteron has only one bound 
state, such labelling has been omitted. 

Hineadly, lt 45fcan@helpatomciarafyatherdescua ptaoni on 


Vij; as an "effective potential": 


Vii j=<on |Ge(E9Ges ‘XK EQGG (Ed gs> 1.48 
=KxF |Gordd Ediixi) 


=<x; | (E-Ho)|x;>. 
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Vo j=<Polg;> Tie 49 
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=<poltEehow xe 


=O on the energy shel] .t 


e. Analytic Structure of the Equations 

The first step toward solving integral equations 11.40 
LSwamscudywotethemsinguleriiyestructure) coma csmeconstituent 
parts. Because many methods of solution involve analytic 
extensions of these parts into the complex plane, the 
following; Drier accountyotesingmlaratmestwadsiainotebe 


restricted to the real q axis. 
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In general, the solubility of such an equation depends 
on the behaviour of the kernel Vr (as discussed in section 
Bereranc Of the solutivon \ rtseli. Of “themtwo contributions 
POmEperhernely, thats duesto, the baré” “two-body t-mattis 7 16 
Ehewsimptlest. In pure wwe body scattering, (Eb) has arsimple 
pole|at each energy corresponding to a bound state of the 
pair. In the three-body case, r(£-q2) will therefore have a 
polepon the seal (o axus foreach E-g2 corresponding tora 
bound state energy. For three nucleons this means at most 
one pole. 

Singularities in the "effective potential" V arise from 
both the three-body propagator Go and the form factor 
<G pr =0 p) Bor most ™calctilations, Brayshaw (B8re68) has 
argued that the latter source may be ignored. Since the 
separable potential, defined through the form factors, is 
only approximate, any reasonable calculation must be 
relatively ansensittive towthe exact form of gy, If the 
Singularveies ine were to playa pant, then this 
requirement would not be met. 

When the energy is high enough that the propagation of 
three free particles can occur, the free Green function 
exhibits two pairs of branch points of logarithmic 
shane se with branch cuts between them. Since the 
positions of these singularities depend on both momentum 
Varvablest. while the awntegration 1S pemiormed over won ly 
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pole 
Brome a 


logarithmic 
Singularities 
from V; ; 


Figure si1l:5) the antegrand of 1b-49ehas singulanities 
in both the variable of integration q” and the 
independent variable q. 


Za 


solution of equations 11.40. The locations of all 
Singularities in one kernel (excepting those due to the form 
paGcOCpeahesett lUuStrated inyllgune llwSeacea function of both 
momentum variables. 

Singularities in the amplitudes X;,; themselves have 
been studied in detail by Brayshaw (Br68) who looked at 
analytic continuations into the complex plane and the 
ultimate solution of equations like 11.40. The complicated 
structure of these amplitudes is contributed by both the 
"effective potential" term and the integral term, including 
those Singularities mentioned above which are independent of 
the variable of integration. In particular, Brayshaw has 
found the behaviour of the X;; to be compatible with 


solution by the method of contour deformation. 


f. Methods of Solution 
Following angular momentum decomposition, equations 
11.40 reduce to a set of uncoupled, one-dimensional integral 


equations of the form 
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Solve this type of vequation given the gsingularity, Structure 


outlined in the previous section. 
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The most common approach is to solve I1.50 for 
mtedratzvon along some contour in the complex o” plane. 
Then, from the solution for complex momenta, the physical 
amplitude, for real momenta, can be calculated. The simplest 
contour is defined by rotation of the coordinate plane 
sCaispemeculvalent Lomintegratcon along the lines co. e. “oaths 
contour clearly avoids the "moving singularities" in the 
Kernelewhen they take thestorm indicated wn Ligure 11 [60 For 
some values of gq’, however, the branch cuts assume the more 
complicated form of figure I1.7. In this case, the simple 
rotation must be supplemented by a detour around the branch 
point at A and evaluation of the integrand on the second 
Riemann sheet (where the contour is represented by a broken 
line). Brayshaw has shown that for small enough ¢, the 
method of contour rotation is compatible with the 
Singularities contributed by the X amplitude itself. 

A different contour has been used by Ebenhoh (Eb72) who 
was able to extract the amplitudes for real momenta 
directly. This avoids the problem of integrating over 
Singularities a second time in computing the amplitudes on 
the real axis from those in the complex plane. 

The presence of "moving singularities" in the kernel of 
equation II.50 prevents a straightfoward solution using 
numerical quadrature. While the changing nature of the 
Singularities would require a changing integration mesh, 
solution of the resulting system of linear equations would 


demand that the mesh remain constant throughout. Doleschall 
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Figure 11.6 Invtheir sinplesb? formpethesbranchgeuts 
in the integrand can be avoided be simple contour 
rotation. 
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Figure II.7 For some values of g the basic rotated 
contour must be adjusted to negotiate more 
complicated branch cuts. 
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(Do73) (and also Sohre and Ziegelmann (So71)). has developed 
a modified Gaussian quadrature which gives exact solutions 
for an N-7Tst degree polynomial using N mesh points. It 
should be noted on surveying these two popular techniques 
that analytical integration inevitably requires that the 
potential have a suitable analytical form with a well 
behaved extension into the complex plane. Numerical 
integration, on the other hand, places no such restrictions 
on form allowing a wider and perhaps more realistic choice 
of potentials. 

A number of calculations have been based on the 
Fredholm method of iterating equation 11.50. A problem 
arises from the fact that the multiple scattering series so 
generated does not always converge. Kloet and Tjon (K171) 
enforced convergence by adding a multiplier to the kernel. 
The transition amplitude, which was derived for a local 
potential, was then approximated using the Padé method. 

By splitting the X-amplitudes into singular and 
non-singular parts and iterating once, Stuivenberg (St76) 
emnectediiassoOlutmionsby analytical integration along the real 


axis gon ly. 
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IIIT. IMPULSE APPROXIMATION 


demJustification 

Although the impulse approximation was originally used 
ioe enrinact Learuress Ol sienean-Mmminberactionma Lroman—D 
inelesticescattering® (Cho0jer4 tseiquantiteatavel accuracy has 
Since been surpassed by full-fledged three-body methods. For 
esr samplierty,) however,;s thecimpulsetapproximation is still 
amusetul tool of qualitative analysis, especially in the 
vicinity of quasi-free scattering. We have used it to get a 
Eee], fore thebsensicivitiess otether Din, nn) p  quasisfree 
scattering reaction and to guide us in the most efficient 
use of the lengthy and expensive three-body calculations. 

Kamenacn  aatliys n-n quasi-free scattering can be defined 
by a three-body final state where the proton has no recoil. 
Experimentally, this state is evidenced by a peak in the 
break-up cross section. Except forthe effect of the binding 
energy, then, this is two body n-n scattering with the 
proton aS a spectator and should be well described by the 
impulse approximation. 

As incident energy decreases, the binding energy plays 
a greater role and the reaction adopts a more three-body 
character. Under these circumstances, the QFS peak is less 
pronounced, broadened by the contributions of various 
Moe leriscat bering ier mom Phe ene royworrin sash uly) G2 1h 
MeV (14.3 MeV in the centre-of-mass), is still high enough 


relative to the deuteron binding energy that the impulse 
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approximation should give a reasonable account of the 
kinematries Of this reactiom and an indication of its 
Suitability for measuring the n-n effective range. 
Furthermore, the impulse approximation has the added 
eQvantade=thierhiteanedudecshe implicitly, aularebatiyeunoD 
partial waves in the incident and scattered channels. 
Stuivenberg (St76) has acknowledged the importance Si higher 


Pertialwwavesyin this particular reaction. 


b. Model 

Dhes2mpulsesapproximation modclmbor tiesOrSereactictm1s 
illustrated schematically in figure III.1. The kinematics 
are those of simple n-n elastic scattering modified by the 
absorption of energy in the breakup of the deuteron (see 
appendix A). 


Yamaguchi S-wave separable potentials (Ya54a) 
V(kK,K’ J=g( k)dAg(k’) igdel’.<d 
where 
g( k)=1/( B2+k2) ee 
are used for both n-n and n-p interactions. The momentum 


space proton wave function in the deuteron thus takes the 


form 
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Figure III.1 In the impulse approximation the 
deuteron wave function (at A) furnishes the momentum 
distribution for neutrons scattered by the off-shell 
None interaction™ (at = B94 * The=protoneactstast a 
Spectator. 
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W(p)=N/[(a2t+p2)(B2+p2)] qlee 


Whebe@gN is a normal izataoneconstomts lt "si the peaking. of 

thisewave function at zero momentum that 45 responsi bles tor 

eherpeake ina the! m-Decrosst sectiontat i thesOrstgqeometry.: 
Porpchemseparabletpotenboalgitizag, heltomatrix for nn 


scattering is also separable (see equation I1.16) 


tlk, k’sE)egl(k)r(E glk’). Wied 


For on-shell (k2=k’2=E) S-wave scattering the effective 


range expanded f-matrix is used 


BU K=129 =CKeoto= 1k) |! td leo 


RMT 2/2oP res 1K 


where @ is the scattering length, ris the effective range 
and P is the shape parameter. Equations III.4 and III.5 
imply hate theshalteofi=eshel letlmatrax at Braattigure Tie 


can be approximated by 


Bibeke 2 key =| ot kee oun (np arene7eeP rs ikse Kyo leat 6 


Here the off-shell effect of the potential on the behaviour 
of the t-matrix can clearly be singled out in the factor 
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The complete expression for the cross section is 


seal he] AlCepyer (eee: Dil. 7 


wneresthes last ‘lactor mS the kinematic factor in the limit 
of zero proton recoil (see appendix A). Although 
approximate, it is considered to be in keeping with the 
aeeuracy Of the impulse approximation;yet least in the 
VEcinity of OFS “geometry. 

The effective range parameters used in the above 


Patveractions are: 


hen p= er OM io ee oak aii 
Copanasinglet) )'a=417.0rmn CSpini str plet)ea=-25. aM 
P= 204197 


c. Experiment 

As this study 1s being conducted in Support of an 
experiment in progress at the Nuclear Research Centre, a few 
words about the experimental setup are in order. Note that 
all angles and energies (except relative energies) are 
described in the laboratory frame of reference. 

Incident 21.5+.3 MeV neutrons scatter off a D20 target 
into two banks of detectors arranged in four opposing pairs 
at +42° to the beam axis with 11° azimuthal spacing between 
pairs. The detectors have an angular acceptance of 320 nee 
Gudirec tron anded2 sin these direction = (her latter value can 


be subdivided into two or three equal-sized bins if 
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Figureslil.2 Schematic diagram of @thesD(n, nn) p 
experiment. Detectors are arranged in two banks 
around a 42- cone. Azimuthal Spacingewithin each 
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necessary. The experimental geometry is illustrated 
Scnhemelically,in, figure a il..2, 

A n-n QFS reaction is evidenced by coincident neutrons 
in each detector of a pair. The proton, remaining at rest in 
exact QFS, cannot be detected. Neutron laboratory energies 
from 2 to 19 MeV will be measured with a resolution on the 
order of 1 MeV. Symmetric QFS occurs at a neutron energy of 
9.64 MeV. The goal is a meaSurement of the differential 
cross section at n-n QFS geometry, to an accuracy better 
than 5%. A value for r,, can then be deduced by matching the 
measured cross section to the value predicted by a 


three-body calculation. 


d. Results 

For on-Shell scattering, equation III.6 reduces to the 
Simple effective range expansion for the t-matrix, Since the 
Off shelt function glk ey/oqlkvaas unity .othiset-matrixy being 
defined by measureablet on-shell quantities, should be 
independent of the off-shell behaviour of the potential. 

lem thei absencesofia clear! choige, for) the, nan potential, 
one should therefore seek out kinematic regions where 
off-shell effects are minimized, if a reliable estimate of 
Po use toubermades» Ofi icounsesithe, exactsOrSsqeonetiya is best 
Suited, but the experimental limitations on resolution and 


acceptable count rate demand an investigation of near-QFS 


geometries as well. 
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POmer igure Ti) ieecne off-shel li runction o Ki /ol Kk) 
depends on the incoming virtual relative energy Eloff) and 
Ene seal outgoing energy Elon) (k* 2=E( off) K2=E(on)) one 
could, therefore, regard the difference between E(on) and 
E(off) as an indication of the model dependence of these 
Calculations. Fom BlOff ictoser to Elon), the reaceion, has 
more of a two-body character and should yield a more 
Belpoblesvalue for ira. Proton Lecol leenengy menoOwaver at cea 
further test\ol two-body|character*and must be mona tored as 
well, lt istalso aimeasSure of the validity of this 
approximation. 

In the experiment, neutron coincidences are possible in 
eetectors weth 169°, (T58@.and 147° azimuthalspaganga(t3 0). 
As the separation decreases from 180° the calculations 
become more model dependent, but not greatly so (figure 
III.3). At 180°’ the cross. section is constant over the 
azimuthal range of the detectors (3° each). For coincidences 
in other than opposed detectors, the proton can acquire 
recoil energies comparable to those of the outgoing neutrons 
although this problem is alleviated by using energy windows 
Lore the neutrons. Figure III.4 shows the magnitude of proton 
recoil energies encountered, for fixed neutron energies. 

In the case of fixed neutron angles and varying 
energies (tigure ILiwoga mode dependence increases away 
from the OFS value of 9.64MeV while the constant on-shell 
energy implies that the cross section will only be sensitive 


forthe magnitude of the on-shell amplitude but not topits 
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k-dependence. Over the range of energy resolution (1MeV) 
model dependence is constant but the peaking in cross 
S€Ct1ON 18 Still discernable. 

In contrast to the above results. the two-body nature 
of the reaction seems to persist through variations in the 
detector elevation angle 6 relative to the beam axis. For 
symmetric opposed neutron flight paths both on- and 
off-shell energies increase with angle (figures,II1.6), ‘the 
nee ekiect being that the intluence of off-Sshel!l potenti a. 
Denaviourn|s constant (figure Tide. 7) This, and the 
generally small proton energies encountered (figure III.8), 
indicate that the 12° angular acceptance of the detectors 
should not pose a serious problem. 

The cross section, however, iS quite sensitive to 6 and 
peaks at an angle of 39°, somewhat smaller than the 41.7° at 
which the kinematically defined QFS occurs. This is clearly 
because of the on-shell t-matrix which increases as @ (and 
hence E(on)) decreases, displacing the peak from that due to 
the deuteron wave function alone. The fact that the cross 
section changes by a factor of two over the angular range of 
the detector will have to be accounted for either through an 
averageing process or, perhaps better, through the binning 
mentioned in the previous section. 

ThessUuppresolonmert Om -snell eErlec Shoversamrangesot 
Symmetric, non-QFS angles, coupled with an almost linear 
relationship between E(on) and 6, suggests a variation in 
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expanded on-shell t-matrix (equation III.5). Figure III.9 
demonstrates that this shape dependence is too small to be 
of any use. 

The relationship between o and r,, is roughly linear 
with a slope of -1 over the conceivable range of r,, (figure 
111.10). From this, a determination of r,, to within a few 
percent seems at least feasible. 

AS a guide to identifying possible sources of error in 
the experiment, equation III.7 can be inverted to yield ras 


a function of o, a, and P, whence 
(Ar)2=[(dr/doa )Ao ] 2+[ (ar/da) dal 2+[ (ar/dQ) AQ] 2 iG 


In the last term, the shape parameter P has been 
incorporated into Q=Pr3. Comparison of common potentials, 
both theoretical and empirical, suggests that there is some 
correlation between P and pr which would permit treating Q as 
an independent parameter. 

For reasonable choices of the above uncertainties, the 
normalization of the cross section arises as the largest 
Source of error for all possible relative outgoing neutron 
energies (see figure I11.11). Even at the reasonable level 
of 5% SOREN OAR ONE the cross section still 
contributes an error inp that is as large as some recently 
reported values (table 1.2). One must therefore hope either 
for a significant improvement in the measurement of detector 


efficiency or for a noticeable sensitivity of the shape of 
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the cross section tor. 

By way of comment on the liberal use of the effective 
rangemexpanded f-matrix in this part, 1£ 15 instructive to 
consider the expansion in terms of the dimensionless 


parameter rk: 


kcoté=(1/r)(-r/at(rk)2/2-P(rk)44+...). III.9 


For reasonably behaved coefficients, rmk<i1 should guarantee 
convergence of this series, whereas the reaction at hand 
yields rk=1.2. While this does not preclude convergence, it 
is true that the third term contributes as much as 7% in 


this particular case. 
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IV, THREE-BODY CALCULATIONS 


a. Model 


Three-body calculations were carried out with the 
Sompurceurcode OfeP.-Déleschalin(Do/3).eOur® versionaer the 
the program package (which we call DOCODE) is written in 
IBM-compatible FORTRAN and runs on the University of 
Alberta's AMDAHL 470/V6 computer. 

The program solves the AGS equations (11.30) for 
SepalebteyNeNrinterace1cnsae(liesos ReThenreducedsintegral 
equations (11.50) are integrated numerically by Gaussian 
quadrature. Those parts of the kernel which contain only 
Simple poles (from 7) are dealt with in a straightforward 
principal value integral. For the remaining more singular 
parts (containing logarithmic singularities and cuts) a 
Special quadrature 1s constructed based on the mesh used for 
the simpler quadrature. The method relies on some smoothness 
Uwe oplittangesof-the kernel antes’ singuiaripend Ezequilar® 
parts is done in such a way that the special quadrature is 
independent of energy and two-body interactions. 

Integration yields a system of linear algebraic 
equations which are solved by the Padé approximant method 
(Ba71) using MacDonald's e algorithm (Md64). The total 
interaction is constructed as a sum of partial wave 
three-body t-matrices up to 19/2* total angular momentum and 
parity. To speed up calculations, orbital angular momentum 


which would normally include states up to £=11 is truncated 
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atecolewith no effect on results reported @netheloriginal 
calculations (Do73). Our calculations show a 2% difference 
between n-n quaSi-free scattering cross sections with 
Gruneation at =7tand) ’=8 airor partrvaliGwaves beyond 7/27, 
quadrature is replaced By Sloan's unitary first-order 
approximation (S169). This is a’type of first-order 
approximation to the multiple scattering series with 
Unotane ye DULIt an. Our calculations smamechercjuasi—tree 
Scautering region mndtcavel an effect ome, in the ecombined 
Oeergana 11/2* partwalewaves and 04% amethne combined 811/2° 
and #i3/2°> partia iwavesi-Sdtle to this Ssubstrievtion. 

Except for the numerical treatment, the distinguishing 
feature of any three-body calculation is the type of 
two-body interactions employed (see table IV.1). For most 
practical purposes, the mathematical complexity of the 
problem demands the use of separable potentials. Limits on 
sheer computing power place further restrictions on the rank 
of potentials, precluding for example an accurate series 
approximation to a, provem Local spotential; or even tne 
inclusion of an OPEP tail. 

Phitial calculations with DOCODE®reported) by Doleschall 
used rank one Yamaguchi (separable) potentials in S-, P- and 
D-wave interactions. Later the program was modified to 
accept generalized Yamaguchi potentials which were designed 
to improve the description of the mixed-state deuteron and 
to fit the experimental phase shifts in all partial waves. 


In particular, the polynomial form of the numerator enables 
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reference reaction 
(Aa65,66) break-up, elastic 
(Ph66) break-up, elastic 
bound state 
(M169) bound state 
slo) elastic 
Geagaly) break-up (FSI) 
(Ki > break-up 
(E572) break-up 
CHET 23) bound state 
(Pi72) elastic 
(Doves (82) break-up 


CD73 ) break-up (OFS) 
(Ha73) bound state 
(Sk75) elastic 
(Be76) elastic 
(St76, 78) break-up 
(Bn77) break-up 
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rank partial 


waves 


Table IV.1 Survey of three-nucleon calculations. n-D 


system with separable potentials unless otherwise 
noted. T indicates tensor force (°S,-*D;). 
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reproduction of nodes in the ‘Sot, °S, and *Py phase shifts. 
Currently the program accepts both generalized Yamaguchi and 
R-type (recently developed by Doleschall) potentials of any 
rank in any partial wave. Coupling between partial waves, as 
Piet nemcsh.=°D yitensomlinteractrome lisifta Vocravanilable. 

The tgeneralized Wamaguchi ‘potent fal & ymethe kK=th 


partial wave, has the form factor 
RG" (pian SP aoy moony [ie eres apes) DVL 


with m no larger than 2 for the potentials presently in use. 
This form will hereafter be designated as Y-type, as a 
reminder that it 1S more general than the traditional 
Yamaguchi potential (Ya54a,b) which inspired it. Of course 
in the case of rank one S-wave form factors this distinction 
ismpure: formalityeyeThe iR-type storm vactonme fon which etiere 


is no published description, 1S very similar: 
Gp y= ip ACE B p aS) 2 tam ley Maes Bee | IV.2 


Vabhiesome iearrangement yl wt teanebe tseenechat (ihe Reform 
farLormns eceatlyearmreduced orm of thewy, withe)>-0Fand 
Be=6. Thewsmaller number of parameters int this storm factor 
(tnvougm its restricted ttorm) “allows @itsiipractical mse? for 


Namen values ot: mi(up to “Si ntour calculations). 


+Standard spectroscopic notation (***'£;) is used. 
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Note the distinction between the number of terms in a 
form factor (m) and the rank of the potential (number of 
Channel svor, in mostecases}« forme factors, ! indexed bysy and 2 
in equation 11.35). This terminology will be used 
throughoutst Furthermore; Ptosstreamline notationenthernrank of 
an) S-waver potential may be andicated by a prefixson the type 
CeSvionallony Wass imino one etce 

Both types of potentials are used in rank one and two 
Singlet S-wave forms, this being the most important partial 
wave. Triplet S-wave potentials (n-p only) are rank one 
except for a rank four tensor interaction (abbreviated 4T) 
used with the 2R singlet interaction. P-wave interactions 
abemexclisavely Yrtyperandm@anchudesthei dP S-ebye Beeeandase, 
partial waves, the latter represented by a rank two 
potential. D-wave interactions are rank one R-type and 
include the ‘Djfand °D2, partial waves. The °D, interaction 
is only present when tensor forces, which are R-type, are 
included. All interactions allowed by the generalized Pauli 
exclusion principle (4+s+t odd) have been used. Traditional 
Yamaguchi and Y-type tensor interactions have been used by 
Doleschall (Do74) but have proven unable to fit even the 
two-body experimental data. 

Parameters for all interactions were obtained from a 
weighted x2 fit to accepted low energy nucleon-nucleon data 
(Mg69) including singlet and triplet S-wave scattering 
length and effective range, Pa, Qu, Au/A, and p-p phase 


shifts up to 460MeV. The latter were used to construct both 
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Mb and iis inkeracthons aby mwereasungithesass1 oned 
experimental uncertainties below 10MeV, thereby diminishing 
whe Ewemgqiessyor cuhe Biowempenergwiphasessiafts. In this 
respect the calculations are charge independent. In fact, 
wherewer both n-m and n-p interactions are present in the 
higher partial waves, they are identical. But since the more 
important S-wave effective range parameters were allowed to 
differ between n-n and n-p interactions the calculations may 
be regarded as charge dependent. 

Interactions obtained from the same two-body input data 
are nominally considered on-shell equivalent although there 
may be on-shell differences due to differences in overall 
form. Effective range parameters used for S-wave potentials 
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all partial wave potentials with the number of parameters 
available for fitting in each form factor. The actual values 


Ob sperameters' can be found im appendix C. 


b. Results 

mie. calculated f-=n OFSECcroSss Sectirons@ ter the varicus 
two-body potentials are shown in table IV.4-6. In the 
abbreviated notation for the potentials the Foner term is a 
M-h Interaction, the second}is n-p. n-p cross Sections are 
included to explore the results of von Witsch et al. (Wi80). 
In addition to the overall differences in form between Y and 
R potentials we had available to us the choice of rank (of 
Singlet S-wave potential), form) of triplet)interaction, 
presence of P and D-waves and value of the n-n effective 
range. Except for the last quantity, the effects of these 
are treated as independent. 

The completet 2R interaction is the most sophisticated 
we have used both in terms of complexity and fit to the 
experimental two-body data. For this reason it is considered 
to be the most reliable. Assuming this, the most striking 
Peaturer of tablevlV 14) 16 sthes bal LULeEmO met hes I Re Interaction 
Pomproduce the correct Magnitude of crossmsection | Thiceis 
believed to be due not to some peculiaritysin thewik 
interaction but rather to the more general fact that 
Symmetric QFS is indeed off-shell independent, but ond vaert 


rank two (or higher) interactions are used. By the same 
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Table IV.3 Number of terms and individual parameters 
avaibable for fitting. Refer Comequations: 1Vieand 
[Veco Ines consuructvon of «tensor 
shown in appendix C. 
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ie Src oe 
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token, the remarkable accuracy of the 1Y potential is 
regarded as accidental although not entirely unexpected, 
given its unusually good agreement (for a potential with 
only two free parameters) with higher energy phase shifts 
(St76). Ignoring the 1R calculations, variations in the 
cross section predicted by pure S-wave calculations are less 
Revet ries wee 

Pesomewhat stronger result is the lack of any 
Spprecciablesisensitivity to thempresencesoffa tensor force, 
exmmbicedrby! both Ri andy 2Ricallcullat ionsel Theil facteathatea 
rank > one tensor force’ was used in the former case while a 
rank four tensor force was used in the latter case 
reinforces the generality of this result. 

The effect of adding D-wave terms is roughly double 
thateore P-wave terms but. oppositeranhe dairectione @Thisikaind of 
interference between odd and even partial waves has been 
observed in other kinematic regions of the three-nucleon 
system (Ku81, Ke78). The combined effect of adding P- and 
D-wave interactions is on the order of 1% (decrease in cross 
section). 

Phesonesteature worwhich ell calcoudarions were 
Sensitauve , and: with’ considerable consistency, was the value 
Dre ien For! avdecrease of 8.1%) an mtn) SYMMELTIG IORPS@cross 
Sections increased by 8.3-9.4%. In determining a value for 
r,, from a measured cross section, therefore, percentage 
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A slight inconsistency will be noted in the use of a 1Y 
n-p potential in the nominally 2Y calculation with 
Preae-O1/M, This was done £o simplify (i.e, cheapen) 
ealeculealions, assuming that this) potential plays) little ware 
Unen=n OFS. Our results show that its effect even on n-p OFS 
1S minimal. Given this, and the agreement between 1¥Y and 2Y 
Cadeulataons anyway, the substitution shouldehot prejudice 
results. 

TiemecrOSse SeCLION Or Symmerric No OFSers ) in tcc ce 
insensitive to everything except the triplet interaction 
which all but uniquely determines the magnitude. 
Consequently ithe ratio of n-petoin-n OFS cress Sections, 
used by von Witsch et al., behaves very much like the n-n 
OEGmcroecs Seeuion alone, For thes8.1% changesin , this 
baedo fella by 8.0-9.1%. It is therefore no more sensitive to 
Pheevalue, off7ree than is the n=-neOrs crossMsection, nores 
its sensitivity any more consistent from one calculation to 
the next. 

The geometry of the von Witsch experiment was actually 
Slightiiveasymmetricyeclose, to thatpotetables1Vl5.eAllaof the 
remarks concerning symmetric OFS, however, can be applied to 
this case as well. This includes the theoretical] result that 
fheatationct 1-p to n-n OFS cross Sectionseis no more 
Peliable an indicator of P,, than the m-m OFS cross iseciicn 
giselt. 

For radical departures from symmetry, as in table IV.6, 


dittferences due to two-body input are much clearer. The 
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Table IV.5 Centre-of-mass cross) sections: for n-m and 


n-p n=pynen 
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Table IV.6 Centre-of-mass cross sections for n=n and 


tii=ye) n-p/ non 
- foe 4.336 
748 Bir pe hone 
belay a a 339 
. 7490 4.049 
27741 4,594 
Peeve 4,476 
. 7648 4S: a Sa 
7Os 65 4,745 
. 7646 4.176 
s hbT.S 4.266 
Ashes 3.654 
.O263 32055 


feo ORS) ate neutron polam angles! OMNlG.2 @.8oF. he. 


(units of mb/sr2MeV) ¢2-¢,=180°. 


Not ineduded, 


Phase space factor 


64 


ez ia Mp 


2.2% 


seat 88 


rsesve — is2 


i Ase met ss 
tae arTe sae 


Saa a9 303 coat 296% tab Ab—rt a6) 3. NS ee 
i idiece:s: as 


Lo a5 44 he. 
eat Sah ie. +r 
- obey Smet). 


i 


7 


etiect ‘ofethe tensor forcesisesiightly enhanced but 


negligible compared to the increase in cross section 


produced by the addition of P- and D-waves. 
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The properties of the break-up reaction at symmetric 


non-QOFS neutron angles, as predicted by the impulse 


approximation, seemed interesting enough to warrant further 


Study with the full three-body calculations. Furthermore 
with the large angular acceptance of the detectors in the 
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Mable iV.) Laboratory break-up Crossmsec clonal on 
symmetric geometries. (units of mb/sr2Mev) 62=03, 
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Pomceocuscsed further ij echepter V: 

Table 1V.7 lists the predictions of the various models 
FOr the break-up cross section with equal’ outgoing neutron 
angles. The most prominent feature Of the results when 
viewed graphically (figure IV. 1) is that the cross section 
peaks around 40.5% (corresponding t0 a proton recoidsenergy 
CheikevV). This 16 not as low as the impulse approximation 
Peeadters) but still signiticantly Jower tham tnemOns uanoremc: 
ae) This 1S important im Wiew Of the steepeadrop mnmearoce 
section above 42°, much steeper than that Peecthasse by the 
impulse approximation. 

The sensitivity to various aspects of the two-body 
input is best illustrated by ratios of cross sections 
(igurel mV. 2). Sensitivity to the rank of potentiale 
represented by the ratio of 1Y¥ to 2Y cross sections, is 
fairly constant’ throughout the angular range of interest. 
The effect of D-wave interactions, measured by the ratio of 
complete 2R calculations with and without D-wave components, 
Piewmeasces MOMOLOnically With angle. From the vesults ac 
Give ercecan pe acsumed that the Carecrumot Pawave 
interactions also increases with angle although smaller and 
Opposite: in wign to the D. Discrepancy betweens ihe 
predietions) of the 2Y and complete! 2Recalculations, 
interpreted as the effect of the form (ore Oru sSnedimecontelon 
of two-body potential used, decreases monotonically with 
angles Overall, there Seems tO. be no Mpreterred ~anglesat 


which undesireable sensitivity to two-body input is 
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minimized. 

Sensitivity to fay secalculated) with the iy potential, 
1S Significantly greater than the above sensitivities and 
rises steadily with angle (figure IV.3). For the 8.1% change 
in ra, the change in cross section varies between 6% and 12% 
over the range of one detector. 

Finally, predictions of FSI calculations with the 2R 
interactions, detailed in our published results (Do82), 
illuminate the different character of this region and its 
UniSuLtapwlity fom a measurement Of 7, ,.s5OEN Nano and sn spears! 
cross sections show significant sensitivitiesf to the 
presence of P and D-wave interactions but not to the tensor 


force. 
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V. CONCLUSIONS 
Ideally the presence of experimental and theoretical errors 
m@che determinakvem of some tquant: ty ishouldiberrefleccedein 
aut Mia menumeriveal Buncentaintyern tehavetqtant uuy dancethe 
experwmentesnould be cameteasout imecuchmaaway mas arc 
minimize this uncertainty. While this should be a guiding 
precept in the present experiment, we see in trying to 
extraceEthis =tworbody «parameter sirom its iciicak lof Sthree=body 
MoLceraahiens lthaaithets wtuaettomens nohisome lear ecuruelne 
experiment, from design through to analysis, must be 
designed not only around known theory but around sources of 


theoretical uncertainty as well. 


a. Theoretical 

Except for a few instances where quantitative evidence 
was available (see section IV.A), we have assumed there to 
bewmne tenrors in Stheanumerical methods, approx ymations Gor 
other aspects of the actual three-body calculations: While 
Ei? Saiia ys inch be .eaue srand! Giscrepanereseoetween isupposedly 
Simulbarecalculations are troutinely oeportedaAm athe 
beberabure, (its assessment) is "certainly beyondeuhe scope of 
thLGework or mikewnserttalt hough@this tpornesiSstseldom 
disputed, $b82) the entire formulation of three-body 
interactions is assumed to be correct (i.e. capable in 
principle of reproducing the experimental data). 
Consequently the greatest identifyable source of theoretical 


Uncertaintyerests iwi thethe choice jot el woebedy Sinpureto tiie 
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sensitivity of the break-up cross section to higher: (P 
and D) partial waves in the two-body interaction is smallest 
SemSyiMetmiceOFSmbuts sti le) cmgnastiicantwatmshes leveleot 
accULreacyurequiredtint thempresenticases Mheiur) unelusion, 
therefore, must be viewed as essential, although no error 
duce .co tbheareform) etcehhist considered. 

ine elect offmhestensore forces istisama Vardy “smala@ but 
neticeabl exeitsiinclusiion ,-\nowever ,@nntroducessiansource of 
uncer tainiya dues nots jusit vo Fitsboti-eshellh behay voursibutiitora 
cCertaune ignorance ofentsscomrectson: sheild behayiourledin, the 
impulse approximation, the break-up cross section is 
directly proportional to the square of the S-wave component 
of the deuteron wave function. This quantity in turn depends 
on (though not necessarily in proportion to) the overall 
normalization of that component of the wave function. 
Comparison of the pure S-wave triplet potentials (P,4=0) used 
in the three-body calculations with the Reid soft-core 
potential (Py=6.5%) shows a difference of 1.8% in peak 
values of the square of the wave function. Aside from the 
obvious difference between tensor and pure S-wave 
interactions, thes actualvamounte oleD stare! abtnibutablel imo 
the deuteron wave function is still an open question. It is 
currently placed at between 4% and 7%, from which we could 
placesaneupper Mimi tliofedimeon the erroreiromythis scurces | 


Tne oObvi Ouse teste or vhs iSeaitorperntorniiwo calculations 
differing only in their values of the D-state probability. 
Hoechst! quantity wsennputidirectly tos themiameinge routine, 
the job would not be difficult (only expensive). 
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PoechapterULVcrt twaststatred tchati-synmet wicmOnoeis 
insensitive to the form of the two-body potential used as 
fongras@uank StwotorGhigher epetentpallic@a reused feThe mairse 
observation regarding this statement is that a survey of two 
Hathenisimi lar tpotentialsteoculd hardly betca Wied 
statistically significant: Nevertheless we have’ some faith 
UreunemconplereS2RGinteraction parelivyebecausemthe 145 
(R-type) tensor force succeeds where others have failed 
(Do74) ta Greproduc ing tthe stwo-bodyedata Piminehermoremit 
seems natural to use the 2R potential with this tensor force 
in order that all even partial waves have the same form. 

The second observation is that symmetric OFS is 
probably unconditionally off-shell independent but that 
reasonable separation of on- and off-shell behaviour can 
only be achieved with rank two or higher potentials. 
Evidence for this 1S seen in appendix C where parameters in 
the first form factor of rank two S-wave potentials are very 
close to those of the rank one potentials of the same type, 
whereas parameters in the second form factor are vastly 
different. In effect the on- and off-shell properties are 
eoneentrated in ‘thesfinstvand second#@fiorm Pacvons 
wespectively. Likewlse tohemvensorepotentralans consumucced, 
according to the prescription of Ghirardi and Rimini (Ghé64), 
by embodying the deuteron properties in the first pair (S- 
anda D-wave)) Of thorm= factors: 

This places) the®assertion of vonewirtsch tetital Me (Wie0) — 


that differences among potential forms can lead to 
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detberencessoi@gupstoel 0 mins vred icbedrcrossEseotionsce- inea 
new context. The calculations cited are based on rank one 
potentials forawhich thelideayofitphase equivalencevhas 
runmtedumeaning.Rourtofrthersix tpotentialssof LBruinsmaset 
ase Bn akues) pliorBinstance, Bueue Bitted onhybtombhe 
scattering length and effective range. Furthermore the 
greatest disagreement with the Yamaguchi potential was 
ebserved wothethen!  lineary andticuadratna"’ fformetactous:. 
rather arbitrary functions which would provide no guarantee 
of reasonable on-shell behaviour. 

The question of rank, or more generally total number of 
parameters available for fitting, 1S particularly important 
in the present case where it is wished to vary Pan without 
violating unitarity or known on-shell behaviour. The 1Y 
potential with only two free parameters is capable of 
reproducing only the scattering length and effective range 
exactly. If one is changed while the other is held constant 
the remaining two-body behaviour varies without constraint. 
If more parameters were present to “take up the slack”, r,, 
could be changed without such undesireable consequences. And 
if the off-shell behaviour is isolated in a separate form 
factoreurtltoo cameremain junatfectedsbyechanges ingen, 

One source of error not considered that is well within 
Gunegraspestems from the usevot p-pephase shifts to generate 
fieneandeicpwinteractionss alt is iwoutheconsider ing sabecause 
Phemveryephnase shifts that arersof interest, those at low 


energy, are the ones most affected by the Coulomb repulsion. 
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Wey dosAnoteassrqneansernonetonehis sourcerbut novesthatephe 
p-p phase shifts could be corrected by established methods 
(D181). Assuming the weaker hypothesis of charge symmetry 
paps ewoulldsproducesbettestnen Imnteérachionstabeleas bretheea 
PatOreeeoUat TON OL en=ne@ancdener! phascesit Gesecnoulauner 
prejudice the ultimate test of charge independence available 
from the experimental value of r,,. 

The above constitutes a prescription for accurate 
paBeee-bodylcalculattonstalong wrth aateelnieerikeitimies of 
that accuracy. There is no such discussion relevant to the 
impulse approximation as this method is inherently inexact 


and has been used for qualitative analysis only. 


b. Experimental 

The optimum placement of detectors for a determination 
Gnas peirom the OrSScrosstsectioneys guidedtiby a mumber tor 
experimental and theoretical factors. Assuming fixed 
azimuthal spacing, protonPrecotlaénergqySisaminimizedeby 
centreing detectors at 6=41.7°. Moving outward from this 
angle incréeasesttheoreticdal sensitivityvtio tr; and decreases 
uncentainteyndue’ rorthentioum Ofetwo body g@nuéeracumonmused 
Cross section, however, drops off rather quickly beyond 42° 
vyoetamnigea redicédtcounteratewand lessecertametye@abeut the 
theoretical shape of the curve. Detectors with +6° angular 
acceptance would incorporate a large chunk of this region if 
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Asideciromeconsideratwonseotecounterates Pargewchanges 
Inec ross) sect ionlas@amtinaticnboteqeometaymaresimportant in 
Eermss Of BpRopersnarmalizationbetadaragato compare 
experamental potntsrimeastred witheftnaiteeresetutionse toma 
miecnetieal different talecnoss®seatienecurverstheatatter 
must be integrated over the appropriate regions of phase 
Space. This demands an accurate knowledge of the theoretical 
cross section as a function of several parameters. When the 
execs secriontesichanginge rapidly, as withevariatiensei 4 
the reliability of such integrals suffers. For this reason, 
notwithstanding other theoretical arguments to the contrary, 
it may be better to centre the detectors at some angle less 
(by a few degrees) than 41.7°. This was done in a recent 
experiment (Gu80) although no explanation was given. 

Experimentally the "differental" cross section can be 
Obeainedibysdiveding thestotaltcross sectionwionee small 
region of phase space by the volume of that region, 
effectively identifying the average cross section over the 
Beqionawitheuheld@iterentab crossmsce tionGattehe centroid: 
Forma Grossesection linéar in all variables this procedure 
introduces no error. In the vicinity of a maximum, however, 
Greanyeueg nonrinsivhichmthercrossaseciionelseconcave 
downwards, such substitutionswilleresult in’ systematic 
underestimation of the differential cross section. Using the 
smallest possible region of phase space (AF2=1MeV, 
AGOMoRG s=omenAG p=A05=4))) , Pehis ecileceSicmeserimat eateaLaupeto 


2ieinethes present experiment. "As thisehas=been shown to 
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Craneiateidiprectiyeintoga comparable error on Pryzn, we see 
the importance of comparing apples with apples. 

In practice it may only be necessary to integrate over 
the range of symmetric (equal @) geometries. The curve can 
besebtained fromvasfit totchemresultesdof three-body 
CateulattonseaFromyiigunesiiliggtand PhRSSant appearsethat 
changes?®inithedchossisectionsinsother directions aresless 
but a firm statement about this must await further 
three-body: calculations. 

If the detectors are subdivided three ways in the @ 
direction, the possibility arises of observing different 
poinmtston the kinematic locus of OFS. For the correct 
choices of energy bins, the two most asymmetric angle 
combinations will be roughly centred on a point along this 
locus. The remaining symmetric and slightly asymmetric 
combinations are definitely off this locus but may 
nonetheless be compared to the appropriate three-body 
calculations. The large angle symmetric combination will be 
the least useful due to the aforementioned steep drop in 
cross section. 

Events that are symmetric in @ may similarly be used 
although their enhanced three-body character may cause 
greater uncertainty due to dependence on off-shell 
propertresrofl: thestworbodyeintelaction sei heglargertmumber «ot 
permutations available in the next-to-symmetric combination 
(169°) means that count rates at this geometry should be 


comparable to those at symmetric geometry. 
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c. Determining effective range 

Pon Will ultimately be determined using the 
relationship between r,, and cross section predicted by 
Bar eerbodyicalculations. Figure V.1 Shows four points on 
Enver curve, calculated with the iY and 2Y potentials 7: 
Fitted to these, by a multiplicative constant only, 16 the 
CULMe=Dredicted by themimpulise approximations As this 
NEegGleccs mrescattering erfects, at 1s) expected ito yield ca 
Sirohely Greater Sensituyvriy sO! (oO mtO fa weAS can DemSecnm Ene 
Opposite 16 found, indicating that the experimental eross 
section 1S an even better measure of r,, than simple theory 
SCC est ss 

In the most exhaustive theoretical treatment allowed by 
this study, there remain two sources of uncertainty. 
Discrepancies among cross sections predicted by calculations 
Witheaitfering ‘forms of two-body input are less than 1%. 
Since only two forms were considered we increase this error 
estimate to 1.5%. Guided by the impulse approximation, 
uncertainty due to the disputed deuteron D-state probability 
is placed at 1%. Assuming these add in quadrature, a 
Meheoretical error” of less than 27% can be claimed. 

If a further "experimental error" of 2% is present, we 
conclude that Pyz, can be determined from a one point 
measurement of cross section (at symmetric QFS) with an 


Uncercvainty of about 2%.) BY Liltingepoints=icomeseveral 


TEouebUrposes Of disclisSsiom thevdifterencialecross section 
is used. As discussed above, for actual comparisons this 
must be integrated over certain directions in phase space. 
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near-QFS geometries this figure may be reduced even further. 
If on the other hand experimental error is placed. at the 
otherwise reasonable level of 5%, theoretical error analysis 


1S pointless. 
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APPENDIX A KINEMATICS 


n-D break-up reaction, whether described by the 


impulse approximation or three-body methods, is governed by 


overall conservation of momentum and energy. Referring to 


Pigere iil.1, for a neutron ancidentealong the 2 ax1s 


momentum conservation yields three equations 


92=q1COS8@ ,+G2COS@ 2+G3COS8§ 3 A.1 
O=q:iSin®@ ,CcOS¢ 1+q2SIN@ 2COS¢2+G3Si1N93COS¢; 


O=q,SiN@,SiIN¢,+q2SiING2SiNb2+GQ3SiNO3SINO3. 


Energy conservation must include the binding energy of the 


deuteron 


(Ey=2.225MeV) which is absorbed in the break-up: 


Oe ode NO dO Came elm & je 


The factor K2/(2m), which relates the energy of a free 


particle 
has been 
the free 
ane most 


avoiding 


to its wave vector (or momentum expressed in fim) 
omitted. For purposes of these calculations, where 
momentum of a particle is seldom needed, the q's 
conveniently expressed in units of Mev'1’2 thereby 


this factor. 


To use the effective range expanded t-matrix (in the 


impulse approximation), one must compute relative neutron 


energies. The usual definitions of relative momentum x and 


reduced mass u 
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K=(k1-kK,)/2 uw=m,m2/(m,t+ms )=m/2+ A.i3 


yield a relative energy 
E=K2x2/2u=(K2/2m)|k,-k2 | 2/2 A.4 


In our chosen units, therefore, we have for the relative 


outgoing neutron energy (in MeV) 


Blom =Gi-0 Sit272 A.5 


and similarly 
Gey a, Wer oh ey Pa, A.6 


To express these in fm-2 for use with the conventional 
effective range expansion we need the factor K2/m 
=41.47MeVFm2. 

For a fixed incident neutron energy, equations A.1 and 
A.2 constitute a system of four equations in nine unknowns; 
hence the expression of the differential cross section as a 
function of five variables, usually the two outgoing neutron 
solid ‘angles «and “one neutron ‘energy. 'Choosing, ‘ass done in 
the experiment, the outgoing neutron angles leaves but one 
free variable which is often chosen to be the are length 


along the curve of possible neutron energy pairs (the 


+We assume a uniform nucleon mass of 939MeV. 


- a 
ace ‘eS 


aw £ 
{ 7 
- 
ra > :© H ‘ = = 
cee =' ims 


7 , 7 — ay 
Zc:¢ ‘noatesu inabbord beagls 


m-@R (0LI>RS Ses 3- (ai raagekatt Pana TS i tee =29ne @ 
qaCi2uen Patogsun cowl. alz Sq 7 eslaginey ays 
- ni gact 21 en VEN Crs Ee cepthesta aia 


1s 


» edb ce ak esipwe te 
Fs : 
Asya ais a = ra 


: 
Pr’ 


89 


kinematic curve or locus). Differential cross section is 
frequentiy described in terms of Ehis eneray variable arc 
length projection). Selecting a point on the kinematic 
curve, by measuring neutron energies, leaves no free 
variables and hence a results in a kinematically complete 
experiment. Choosing, instead, the conditions of OFS yields 
a different™locusewhich defines pairs of outgoing neutron 
polar angles and energies. As to the azimuthal angles, the 
requirement of zero proton recoil fixes their difference at 
180° (coplanar geometry) and the azimuthal symmetry of the 
problem makes their sum irrelevant. 

The kinematic factor used in the impulse approximation 
(equation III.7) comes from an integration over the final 
State (outgoing) phase space. Following Taylor (Ta72 p.348) 
and recalling that in our units free momentum is expressed 
as /2mq, we begin with the general result (in the laboratory 


frame) 


a(aG2,ds-D be Gay A.7 


(20 )4(/2m)7/2q4 fdQ2q2 2dq2d23q32dq35(E2-E.-E2-E;3 ee 2 


where D indicates the deuteron target, p indicates the 
proton at rest in the final state (under the assumption of 
ZeLompLrocon recoll), & sand feauerthe ino alerancdelina! 


state energies: 
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EY sB I= 64 B=EYfEe te 3° A.8 
Ror spurposes of tthe! integration the f-matrix is assumed 


constant. Using f=q%, \db=2qdquand: divirdingmby tthe 


appropriate differential phase space volume 


G07 0%208),d0E>= iS 


(25 )4(V2m)7/2qh (qg/ 20 .ga7 2 dbs 8 (ESE .-E2-E3)\t|2. 
The remaining integral serves only to determine the second 


outgoing neutron energy Q2, which is uniquely defined by the 


chosen conditions. Consequently, 


do /d2,dQ23,dE 2=( 20) 4V2(Vm)7(qzq3/qi)|t|2. Me eb 


APPENDIX B DOCODE MANUAL 


a. General Features 
The program analyzes a given three-body reaction using 


the AGS form of the Faddeev equations. The two-body partial 


wave potentials must be separable 


aCe VER. SenGe hia, sera e) | B.1 


Wherelsumma thrones Sompl todther rankeo i aenes potentualamhe form 
factor, g, has either the generalized Yamaguchi (or Y-type) 


form 


Gap J=| pie Teoan Oedn 7 | li od it Bie El Bae 


or! thes so-called R-type form 
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where k is the angular momentum. Thee latter’ typevhas been 
recently suggested by Doleschall as it orievidesias cocd Hit 
Pomexperimentaledatal includingt nodest in WenermSio,at oneencle Po 
phasemshulis 2 Poxri thi siiorm, eDolescha Ul siandsam=5 to be 
adequate in the three-nucleon calculation. Given the quantum 
numbers of the interacting particles and the geometry of the 
reaction, the program can produce cross sections and 


polarizations, for elastic scatteringyandethree particte 
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breakup. The program can also be used just to generate 
kinematics. Results are available in both the laboratory and 
cam.) frames, 

Ase OMpletea sourcenlistingr wsecontaimedaum theetale 
DOCODE. An index to the programs and subprograms in this 
Listing@ iss providedvateithe ends cfithisimanualt Mainuprogram 
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b. Outline 

The calculation is distributed among various main 
programs which are run in the sequence indicated in figure 
B.1. Global communications iS carried out via sequential 
files assigned to I/O data sets 1 and 2. These are referred 
HoOkasitatless# andi 26 Inladditionyvecertainmotheruriless 
described below, must be preserved between steps of the 
Calculacd ont 

The GAUS program reads Gaussian mesh points from data 
Set 2 (data file GAQUAD, included in the program package) 
and prepares file 1 for use by the Gaussian quadrature 
routine. Since this task is independent of the specific 
problem being programmed, ataneedibe pungoniy once, forithe 
implementation of this calculation on a particular computer. 
GAUS must be compiled under *FORTX. 

The normal starting point of the calculation is QUAD. 
It separates the kernels into singular and non-singular 
parts and performs a special Gaussian quadrature for the 


singular part. Mesh points for this quadrature are created 
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Nerewand written into file 1. 

KERE deals with two-body interactions. Using the 
Parameters Of “the two-body forum wtactors it catculares 
two-body propagators and t-matrices. These are placed in 
Rulon OMerUntIrTe tuse. 

The pair of programs, FAMA and PADE, must be run once 
POmecach Combination of total Spee backs angular momentum 
and parity. FAMA calculates the kernel of the (reduced) 
integral equation using the special quadrature performed by 
QUAD and the mesh points for the regular quadrature, also 
supplied by QUAD. PADE solves the resulting set of algebraic 
equations using the Padé approximant technique. The 
three-body T-matrices thus obtained are placed in file 2 for 
unlimited use by the following two programs. PADD is a 
double precision version of PADE. 

ELAS computes elastic scattering cross sections and 
polarizations with or without the inclusion of Rutherford 
scattering. BRUP computes breakup cross sections and 
Polanizations forga given) geometry) OgS: andlPois geometries 
can be specifically selected. 

Panally, there 1S Jan auxiliary sprog rally le winch 
takes experimental two-body data and produces, by a simple 
Search technique, the parameters of* the GB-BOhd 
interactions in the form required by KERE. As this program 
does not communicate directly with KERE via the global files 
(eandeoeethe results OLvd Elli erunsmuse be transferred 


manually to KERE. FITI does, however, retain results on data 
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AUS————~=<—_——-2 (GAQUAD) 


1 (QUAD) 


(for each Jz) 


PADE 


(-RES-) 2 


Figure B.1 DOCODE consists of several main programs 
which are run in the sequence indicated by the heavy 
lines. Programs communicate with each other via data 
sets 1 and 2. Other data sets are used for program 
control (5), diagnostics/results (6) and storing 
intermediate results. 


a5, 


set 2 for its own use and can be run repeatedly until the 
desired covergence in the search routine is achieved. 

Phegbulkidoiwishe execution of each main program 16 
carried out in a subroutine which has the name of the main 
programmprefixed by iDQ@e(eiq. eDOGAUS scDOOQUAD™ ete. )) Thegmain 
programs are used to read the TEST array and to initialize 
various parameters. Since some of these will change 
aceording to the’ problem being considered or (in thescaser or 
FAMA and PADE) the particular step in the calculation, 
editing and recompilation of main programs will frequently 
pemnecessary. (Thisisisthe@sgreasomaion divorcing the main body 
of the calculation from the main programs. 

Many of the variables assigned in the main programs are 
checked, during the calculation, for consistency among 
themselves and with the input data. These variables are 
listed here along with the explanation of the input. If an 
error is found the program signals (via the output file), 
providing the required value of the variable, and stops. 
Thus pawhenarn acubteytasquickerun of the program ns 
suffichent to establish the valuevof any such» parameter. 

Subroutine HIBA is the universal error flagging 
routines It DEINtS al Error) code snumber soneuhe OCulpUteL ule 
anomctops the calculation. [he ehrogscodemlOusel si uCcoOmrec: 
parameter as above is zero. 

The complete program package contains a number of main 
programs and subroutines which are not needed for our 


purposes. Questions about these and about the package in 
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general should be directed to its author: 

Drier. Dolescha il 

EentraltReséarch Institute for Physics 

H-1525 Budapest 

Box 49 

Hungary 
J. svenne=at the University of Manitoba’also has a’ good 
understanding Of* the! program,Walthoucgh hisiversion, 1s) pothas 


recent as the one described here. 


Cree lnput 
Input data are listed here for each program, with brief 


explanations following. Each line of the form 
VARIABLE LIST (FORMAT ) 


represents one line or card of input. The format is given in 
FORTRAN. Repeated lines are indented, with the incremented 
index variable given in FORTRAN DO-loop style. 

The TEST array is used to choose various program 
options! Gbach#elementamustibe given the value 0 or 1. 

The FORTRAN data set reference number for input is 5. 


Citpuewis producedion data sets: 


Main Program PRGAUS 


The only input to PRGUAS is fixed and contained in file 


GAQUAD, which must be assigned to data set 2. 
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Main Program PRQUAD 


TEST (Carer X)) 
(MASS (1) )I=1,3) (6(1PE11. : et 
ei Cl) b= 1-3) CT os 2x), 12) 
UGAARGE Cl) ls 1,0) Wii ere x ei) 
Bepeatmoued => NTE 


_ (TWO(I,J),1=1,DTWQ) VAD Va eakes 
(if TEST1=0) SYM,LLMAX,N1,M2,NG1,NG2,NG3,K,L,TLMM,TC1,TC2 
(OLi2ee i) ns Ct 2 bea el ee Se Oe oe ea 


(if TEST1=1) LLMAX,N1,M2,NG1,NG2,NG3 Cig ate Ra) rl eg) 
PEGE a ONLY LLMAX,N1,N1,NOQ,NS MUST BE READ. 
TEST(2)=1 THE WEIGHTS OF SPECIAL QUADRATURES ARE WRITTEN 


ON THE LINEPRINTER. 
TEST Bs =a THE INTEGRALS ARE WRITTEN OUT. 


NiGetovcalynumberrot two-body channels, taking into account 
the ranks of the separable potentials (number of TWO 
cards or lines) 

TWO two-body quantum numbers 

TWO(1,J) specifies interacting pair 

TWQO(2,J) number of form factor in the separable expansion 
(OP fOr a rank ones interact lonem eee 
forra’ rank N interaction) 

TWwOCesd Mangular momentum of pair (X2) 

TwWO(4,J) parity of pair wavefunction 

TWO(5,d) number of L,S pairs 

TWO(6,J),TWO(7,J)=L,S (X2) 

TWOCE) od), TWOKGY JE Second? LE Shpan ren] [Mut Pnone) 

PE there, arevidentical  partacles¥radentacallinteraccions 

need not be duplicated in TWQ; the program will generate 

these. 

SYM=-1 for identical fermions 

1 for identical bosons 
Ostor dHstince parereles 

CLMAR=max[(L Sf&®pairili+[loottpaitirelative tomthitdeparticie) 

Ni number of Gaussian mesh points for P,D waves 
(6 for three nucleons) 

M2 number of Gaussian mesh points for S waves 

(9 for three nucleons) 
N1SM2 

If numbers are chosen which produce disallowed 

mesh points DSQ is printed and the program stops. 

NG1-3,K,L number of mesh points 

TLMM average of two largest L's 

THs TO20speciiyra transformation on the mesh points to 

smooth out branch points in the T matrix 
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Internally checked variables 

WSD1 

WSD2 

N2 

NTC number of two-body channels (rows in TWQ) 
DTWQ number of columns in TWO (=9) 
TEMAGGavenage ohetwoe@largest Mus 


Main Program PRKERE 


TEST Cin ose) 
G2ENTEST 1=0) (NEE tINER2 COG er) ) 
(PepreST =) (or TEST2=1)) TC3TCe Mm (Rael ele a) oo Ah Hira 
CAROL) tien] REE CoE a 2 ))) 
repeat for I=1,NTC 

NGHPAGI, liird=4 2) (2 (252) ) 

repeat for Jsi ATWwo(s7)) 
(BETAKK, Gm lh) eka 1 AT) ChE 12.67 2x) 
(GAMMA(K,J,1),K=1,N2) (HIS SS Pe) 
FEAMBDAWS 70 0a i= Nah) (S(Bise6, 2x) 
ees ea CEQ, OXNEdS C6 FO RETR Ss” eo Saal ) 
(LIM(3 LT), 1=1,NTG) C6\CE iin ae eS) 
(NE (I ? T= NEC) CUT 2s oe tees) 


GEST G1) =" THE FILE CH2 IS EMPTY AND THE TWO-BODY AND 
ENERGY DEPENDENT PARAMETERS MUST BE READ, 
OTHERWISE THE NFE1,NFE2 PARAMETERS MUST BE 
READ AND THE NEW INFORMATION IS WRITTEN 
AFTER THE NFE1,NFE2 POSITION. 

TESPC2)=4 NEW TWO-BODY AND ENERGY DEPENDENT PARAMETERS, 
OTHERWISE ONLY THE ENERGY DEPENDENT PARAMETERS 
MUST BE READ 

EST (34 THE FORM FACTORS AND THE CFF COEFFICIENT WILL 
BE WRITTEN ON THE LINEPRINTER. 


NFE1,NFE2 not generally used, see TEST(1) 
mGserGe are likesyTci tee 
M 
EBO(J) bound state energy corresponding to interaction J 
(Hine eieota two) 
<0 iP unbound 
The target must be bound: EBO(EINF(1))>0 
CIPA coulomb interaction parameters 
CIPA(1,1)=0 (coulomb option does not work) 
CIPA(2,1)=m where m>0 specifies a generalized Yamaguchi 
form hacvor 
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and m<0 specifies an R-type form factor 
OC pY=| pS MiB p eI THe. Aim MOles pe) Wi 


BETA,GAMMA parameters of the form factors, in the same order 
as TWO 
L=TWO(4+27,1) 
Le4/2+m+ 1 if Sm20 
L=1 if mad 
N2=m+1 if m20 
ane) st m<0 
LAMBDA strength of the potential 
N1=DIM(I) 
EINF energy information 
EINF(1) line in TWO of cae channel 
EINF(2) three-body c.m. energ 
=[m (target) /m(total) |E(1ab) EA bandana) 
EINF(3)=N1 maximum number of special mesh points 
(from QUAD) 
EINF(4)=N2 even number (The principal value integral 
requires an even number of points around the 
Singularity) 
EINF(5)=N3 
EINF(6)= 
EINF(7)=NFP 
LIM(3,1)=0 if third momentum interval extends to infinity 
NF(I)=6 or 9 
number of special mesh points in first interval 
for each two-body channel 
The NF(I) can not all be the same unless QUAD uses 
identical quantities (program bug). 


Internally checked variables 

NTC number of two-body channels (rows in TWQ) 
DTWO number of columns in TWQ (=9) 

NLS number of L-S pairs 


CCN number of coupled channels 

NG 

BD 1 maximum number of 6's in any interaction 
GD1 maximum number of y's in any interaction 
NFP Size of breakup propagator matrix 


Main Program PRFAMA 


TEST (er ( xX) ) 
NG,JT,PAR ey ao, 
(if TEST3=1) (FFIND(1I),1=1,NTC) Ges 2x), 12) 
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SUBMATRICES AND THE INHOMOGENEOUS TERM 


LS PRIN TEDS 


ve THE 


COULOMB DISTORTED SUBMATRICES ARE NOT 


MULTIPLIED BY THE CHANNEL COULOMB PENETRATION 
FACTOR, OTHERWISE THEY ARE MULTIPLIED FROM LEFT. 


)=1 THE 
THE 


INTEGER ARRAY FFIND(NTC) MUST BE READ. 
FFIND(I)=-2 ROWS AND COLUMNS AND THE 


FFIND(I)=-1 ROWS ARE NOT CALCULATED. FROM 
FFIND(J)=-1 COLUMNS ONLY THE ON-SHELL ELEMENTS 
OF /THRIFFIND(1)=2 ROWS ARE CALCULATED, THR iJ 
SUBMATRICES WITH FFIND(I).GE.1 AND FFIND(J).GE.1 
OR WITH FFIND(I)=2 AND FFIND(J)=0 OR WITH THE 
REVERSE ARE CALCULATED. 


y= 

1 TWO 
THE 
THE 

yer THE 

a THE 


ONLY THE LONE =STEPUITERATION 1S CALCULATED. 


PARALLEL CALCULATIONS ARE PERFORMED 
PERS [owt THe SCOULOMEED BS TORTRONG 
SEGCONDeW LEHOUTSTHERSGOULOMB DISTORTION. 
TABULATED  PORMPRAG TORS ARB SUSEDE 
BARLIER INTERRUPTED CALCULATION WILL BE 


CONTINUED (IN THIS CASE THE DATA SETS 3 AND 8 
MUST BE CATALOGED). 


NG number of Gaussian mesh points 

JT total three-body angular momentum (X2) 
PAR total three-body parity 

FFIND(I) weighting of interaction I 


Internally checked variables 


NTC 
DTWOQ 
NLS 


number of 
number of 
number of 
number of 


two-body channels (rows in TWQ) 
columns! TWO (=9) 

LeSPpalns 

coupled channels 


maximum number of $B's in any interaction 
maximum number of y'S in any interaction 


=BINF(3) 
=BINF(4) 
=BINF(5) 
=EINF(6) 


=N1+N2+N3+N4+1 
size of break-up propagator matrix 
put=1 if no break-up calculettons sare to be made 


number of 


dimension 
number of 
channel 

number of 


three-body channels (rows in THQ) 


of kernel es 
three-body channels compatible with incident 


open outgoing channels 
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If necessary, make a short run to establish NCH, 
a second run to establish DM, EON and NA. 


Main Programs PRPADE, PRPADD 


TEST A 2k) 
Peron Gis oe, ngs i 
GVETTEST?2=1) FIDI Gi) 
BeDCa EE LOr W=) /F1 DI 

ATOR) Gear i=) NTC) C171 bo) el) 
TES Ea) THE RESULTS ARE PUNCHED. 
TES 2 =| THE INTEGER FIDI AND THE INTEGER ARRAY 


TOFI(NTC,FIDI) MUST BE READ AND FIDI NUMBER OF 
CALCULATIONS WILL BE PERFORMED WITH 
CARRUN GLV=TOFT Cie Lai NTO a= (ori. . 

PEST C2 )= | THE, RURG SOLUTION WILLABRUCARCULATED: 


PES) = 4 THE OFF-SHELL T-MATRIX ELEMENTS 
FOR THE BREAK-UP ARE CALCULATED. 
Pere Sh) =a TWO PARALLEL CALCULATIONS WILL BE PERFORMED 


Cie BI RSD RW PTE eCOULOMBED STOR TRON, 

THE SSE CONDAWEATHOUT VI COULOMBEDISTRORD ION. 
TEST(6)=1 THE CONVERGENCE WILL BE CHECKED ONLY ONCE, 
TEST(7)=1 THE BREAKUP SIMALTRI CES ARERNOTVERINTED: 


I maximum number of iterations of Padé approximation 
=odd number 

TOL tolerance of approximation 
eonditnon 26 24P(n 0) PU ntl oly IATOn| Pel ne Cries ia 

For three nucleons the approximation should converge in 24 

iterations. 

FIDI number of different calculations to be performed (with 
different potential weightings) 

TORT (iva }=2hUIND (ir) ei ortbhesdtbhecaleculacgon 


Internally echecked variables 

ITM =I 

NTC number of two-body channels (rows in TWQ) 
FIDI see above 


NT 

BP2 

Dix 

N1 . 

EQN number of three-body channels compatible with incident 
channel 

NCH number of three-body channels (rows in THQ) 

NA number of open outgoing channels 

DM dimension of kernel 

BR2 

BR3 
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CEQN 
SDM 
SEQN 
BES 


Main Program PRELAS 


TEST (Gr 12, 28) F) 
NFE1,NFE2 (202 20) 
LINM,LFIM,NRM CBee.) )) 
(ie TES@Ms=i))) 1 (Hie) 

CUUNDIGI)) Tae Toi) GRE) Fi 223) 
GCHA RGE G13 )s l= <3) A Gee ee) -ti2)) 


TRSpTC))=1 THE DIFFERENTIAL CROSS SECTIONS AND THE 
POLARIZATIONS ARE CALCULATED. 

EES C2) =a THE NFE1,NFE2 SET IS CALCULATED, OTHERWISE 
THE LAST SET. 

Teor es = ONLY A CERTAIN SET OF T-MATRICES IS INCLUDED 
(TDIM AND TIND(TDIM) MUST BE GIVEN), OTHERWISE 
ALL DIFFERENT T-MATRICES ARE INCLUDED (THE FIRST 
OF THE T-MATRICES WITH IDENTICAL QUANTUM 
NUMBERS). 

TEST (Ayia THE T-MATRICES ARE NOT PRINTED. 

“EE SHIR Sy) i J NEW CHARGES MUST BE READ AND THE RUTHERFORD 
IS ADDED ACCORDING TO THE NEW CHARGES. 

RST (6a THE ELASTIC OFF-DIAGONAL T-MATRICES ARE 
SYMMETRI ZED. 


NFE1,NFE2=1 usually, see TEST(2) 
LINM2maximum value of initial orbital angular momentum 
of projectile with respect to target 
LFIM2maximum value of final orbital angular momentum 
Of ‘projectrle with tespect) tolitanger 
LRM number of points from zero at which the Rutherford 
cross section is to be excluded 
I=TDIM number of matrices used 
TIND(I) indicies of matrices used (matrices are indexed 
in the order that they are calculated by PADE) 
CHARGE(I) 


Internally checked variables 

NTM 

TDIM number of matrices used 

NTC number of two-body channels (rows in TWQ). 
DTWO number of columns in TWQ (=9) 

NLS Mmombe uot uo. Danas 

CCN number of icoupled ‘channels 

BD 1 maximum number of £~'s in any interaction 
GD1 maximum number of y's in any interaction 
EONM 
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NCHM 
NAM 


Main Program PRBRUP 


TRSE (iC Cle 2x) 
CSPR,POPR COAG er Pas 
NFE1,NFE2 (2.1 122K) ) 
LINM,LFIM, SIG (34 i252 X)\) 
MA,NX,MJP,MLP,MPA C Si ces a) 
Gr iLerest3=i )\aTDIM (13) 
(TON DA ty) el =e ana) (Gl Oley Seren 

(if TEST8=1) repeat: for dJ=1,MA 
CSP RVG1 Sey eels 2) (6 (7 PBI 4a) 2x) 

(if TEST6=0) repeat for J=1,MA 
CANGLES(}ee,1=1.4) C6 GIPE 114925) 
(1£ TEST6=1 and TEST8=0) (PAA(I),1=1,NX) CG GPE dia 255.) 


repeat for J=1,MA 
(ANGLESi(ts J)¢ Tin) Cote raties 2x) ) 


TEST ( 
TEST ( 


THE POLARIZATIONS ARE CALCULATED 

THE NFEIGNFE2 SET IS CALCULATED, OTHERWISE 

THE AST LSET? 

PEST (3H =i ONLY A CERTAIN SET OF T-MATRICES IS INCLUDED 
(TDIM AND TIND(TDIM) MUST BE GIVEN), OTHERWISE 
ALL DIFFERENT T-MATRICES ARE INCLUDED (THE FIRST 
OF THE T-MATRICES WITH IDENTICAL QUANTUM 
NUMBERS). 


Tiel 


NO = 
— ~~ 
est 


REST 4 en THE T-MATRICES ARE NOT PRINTED. 

sy oy) =a THE COULOMB EFFECTS ARE NOT INCLUDED. 

TEST (6) =i THE INPUT AND THE RESULTS ARE GIVEN IN THE 
C.M. SYSTEM, OTHERWISE IN THE LABORATORY SYSTEM. 

REST) =< THE TABULATED FORM FACTORS ARE USED. 

TEST OB) a1 IN CASE OF TEST(6)=1 THE QUASI-FREE SCATTERING 
OF THE PARTICLE MPA IS CALCULATED, OTHERWISE THE 
LABORATORY CROSS SECTIONS. THE ARRAY SFKV(4,MA) 
MUST BE GIVEN AND THE KINEMATICS IS TAKEN 
BETWEEN THESE LIMITS. 

LESH o) =A ONLY THE KINEMATICS IS CALCULATED. 

TEST(10)=1 THE ARCLENGTH PROJECTION, OTHERWISE THE 


BOUIDISTANCE STEP IS USED: 


CSPR number of cross section printouts required 

POPR number of polarization printouts required 

NFE1,NFE2=1 usually, see TEST(2) 

LINM>maximum value of initial orbital angular momentum 
Ofmpsevectsie with respectsto targer 

LFIM2maximum value of final orbital angular momentum 
of projectile with respect to target 

SIG=-1 no Coulomb distortion factors 

MA number of angle pairs 
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NX number of points on the kinematic curve 
=-odd number to put quasifree geometry at centre of range 
2S5NXS80 
MJP,MLP indicies of detected particles 
MPA index of spectator particle 
TDIM number of matrices used 
Tae PREY indicies of matrices used (matricies are indexed in 
the order that they are calculated by PADE) 
SFKV Kinematic limits 
There are three possible kinematic loci (figure B.2). SFKV 
specifies two energy windows in terms of arc length along 
the kinematic curve (MeV). The second window can be ignored 
Dy scnocsSing tan interval. ofa G10 °°), 10%) in) cases “A anceb sor 
dmmrcerval of “avarve) im case’ C.. If Hn doubt..a -orelimimany 
run with a large window can help to establish the type of 
curve and the required intervals. 


Internally checked variables 


NTM 
MA number of angle pairs 
NX number of points on the kinematic curve, see above 


TDIM number of matrices used 

NTC number of two-body channels (rows in TWQ) 
DEWO number of columns’ in TWO (=9) 

NLS number of L,S pairs 

CCN number of coupled channels 


BD1 maximum number of B's in any interaction 
GD1 maximum number of y's in any interaction 
N1 

NFP size of break-up propagator matrix 

NT 

TMF 


Main Programs PRFITI, PRFITE 


TEST 

MASS 

CHARGE 

Poke ACH y TOL 

Bepeateror J=1,NITC 
Grwo Ui d)) 1 = 1 pDTWO)) 

KPEO Glee ht NTC) 

(NOC Ue ui, NTC) 

(ENERGY (1 ),1=1,NT) 

(NOO (1) sna NTC) 

I 

(SOFA (153) ,J=1,NFF) 

repeat for NV=1,NTC 
INP, ITM,NST,J,1CH 
SEPQ(1),SEPQ(2) 
SEPO(3),SEPQ(4) ,SEPQ(5) 
repeat for J=1,2 
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Figure B.2 The energy 
forms. 
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Ee 


Ei 


locus has three possible 


ee 


_ (PHAS(1I,J),1=1,NT) VO(T PE lated = oxy) 
(if TEST3=0 and M=2) (SDM(1),I=1,NT) (6 GIPH Ia. 42x) 
(if TEST3=1 and M=1) (PHERR(I,J),I=1,NT) (6(1PE11.4,2X)) 
(if TEST3=1 and M=2) Bepede MEOr 8) =dia2 

OPHERRGESAI)% l= 17M) GOCTBE 1 PRea2x)) 

(SDM Vl) Sra, NT) Co(PE i 4 2x4) ) 

(MERRGT)) P= 12 NT) C6. GAPE iy aos) ) 

GEVA (iE) Sl =i NIT) CHT Xi larete es) 
CPPVAGH), bai 3) G6 (4Betli s &,2%)) 
CW CIs) elk ae) COCGIPEAM -4 250%) 
(VOP(1I),1=1,NFP) RG 9b .G) 5055) 
GSTP 61) “i= PANED } M6 (IPE ea 52x) ) 
CIM GL ) P= e2aneP) CG (GPE iM s4esex))) 
(FMFP(K,J,NV) ,K=1,NPP) CG UE iat xe) 
(FMF(K,J,NV) ,K=1,NPP) (ou iheed 721) 
(FMFE(K,J,NV) ,K=1,NPP) CG Caiile ae 2m) 
(BETA (hed WNW) Ka tip) (oiChit Ss 26,2) } 
(GAMMA(K,J,NV),K=1,NPM1) CS CR iS are Roa) 
UBAMBDAI ENV te, = tN) A( KAPPA Ci -NV) d= 1,M) CSE So 2k )) 


TEST O14 
TST (2 i= THE BETA(O, 
TEST (2) 

CALCULATED. 


MASS in nucleon 
CHARGE 


Masses 


y= 
! j . 


I=NTC number of two-body channels 
J=DTWO number of columns in TWO (=9) 
K=NT number of experimental energies 


L=NG 
ICH=1 


if input (starting values) 


(from a prévious :FITI -run) 


, 


ae 


THE FINAL PARAMETERS WILL BE PUNCHED. 
5 PSBETA (Be 
THE ERRORS MUST BE READ AND THE USUAL CHI2 IS 


(rows in TWQ) 


bs G@inom dataktser 2 


TOM, fi Bt rng Gtolerance 


TWO two-body quantum numbers, 


(never used) 


see PROQUAD 


EBO(J) bound state energy corresponding to interaction J 
Ghaneld tote TWO) 


<O0eTh unbound 


NOO orthogonality indicies 

ENERGY experimental energies 
NOO(I)=CIPA(2,1)=m for interaction I 
I=NFF number of points (energies) for fitting 


FOFA energies 


INP index of incoming particle . 
ITM maximum total number of steps in fit 
NST maximum number of steps for any one parameter 


wroverrides TES 
ICH as above 


SEPQ(1 
SEPQ(2) effecti 
SEPQ(3 
SEPQ(4 
SEPQ(5) asympto 


TZ 


ve range 


tice DS 


scattering length 


) 
) 
) D-state probability 
) quadrupole moment 

) 


(row 1 of) TWO) 
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PHAS experimental phase shifts 
PHERR errors in experimental phase shifts 
The error is used to determine the weight of each phase 
shift in fitting to the experimental data. In order to fit 
nn and np interactions with pp data, the low energy phase 
Shifts are given artificially large errors. 
SDM mixing parameters 
MERR errors in mixing parameters 
EVA(I)=1 if the data for ENERGY(I) is to be included in the 
chi-squared 
=0 otherwise 

PTVA(I)=0 phase shift turning values (not used) 
Weewvelght, faaborsefore kit Aimg 

WF(1) phase shifts 

WF(2) mixing parameters 

WF(3) D wave 

WF(4) scattering length 

WF(5) effective range 

WF(6) D-state probability 

WF(7) quadrupole moment 

WF(8) asymptotic D/S 
VOP order of parameter variation 
The indexing vof Pparametersh using the motatvon Ds ® APaNny, 
(see equations B.1-3), is as follows: 


| Zz 3) 4 2) 6 i 8 9g 10 a 


Bowe payee ones see a4 A4V2 Boee ee, yi ire r22 


If a number outside the possible range is encountered, 

PEbtanquereturns to the first) parameter, sl ieaezenopic 

encountered the remaining parameters are ignored. 

STP step size for parameter variation 

LIM limits of parameter variation, in pairs (upper lower) 

zero indicates no limit 

FMFP arguments of the form factor 

FMF values of the form factor 

FMFE errors in the form factors 

BETA, GAMMA, LAMBDA starting values of the potential 
Dolaanere ss 

KAPPA multiplies the form factor 

=ienonoma bly 


Internally checked variables ; 
NTC number of two-body channels (rows in TWQ) 
DTWO number of columns in TWQ (=9) 


NLS number of L,S pairs 

CCN number of coupled channels 

BD 1 maximum number of 6's in any interaction 
GD 1 maximum number of y's in any interaction 
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d. Sample Run 


The following control files demonstrate the n-D OFS 
calculation denoted by 21S ,Ry2'"S GR 4Ar4RAseP, Divine 
calculation, the most sophisticated we have presently 
available, uses rank two R-type. singlet Sete Ge potentials 
waithi ‘al imank touricensor tonce @ndmhighertipartial waves (P,D) 
added. The complete interaction, in the order specified in 
Live TWO matrix je could: be deseribedi by: 

aligars BOG pe Oe Len et pe aie 
Np mee San Paleo Poy 8 hav) Geo LUO Ree Doe penetra 

The complete control files are reproduced both to 
illustrate one convenient method for editing and running the 
main programs and to provide the necessary values of 
internal program variables and array dimensions, 

The procedure followed is to copy each main program 
from the master source file of main programs, DOMAST, make 
the necessary changes via the system editor (still in batch 
mode), and compile into a temporary object deck. In this way 
the master file is never actually tampered with. Note that 
the batch output echoes the original tile, not the edited 
version. The subroutines, which never change, Bene: in the 
(don-library) file DOLIBR, insobject torm: 

The three-body T-matrices are calculated up to 19/2* 
three-body angular momentum and parity .eForm this, the: usual 


partial wave limit for three nucleons, several hours of time 
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on the AMDAHL are required. Higher states must be added in 
Domes Gt (22.172) A eG) 2 oy eC arG Intel ray aan 
desired, the cutoff must be reduced by pairs. 

Up to 9/2*, certain parameters are constantly changing 
and JOBFAPA is used to run individual iterations of 
FAMA-PADE. Beyond this limit Sloan's approximation is used 
and internal variables remain constant, allowing the 
remaining iterations to be performed by JOBFAIT. 

The use of GAUS is not shown. AS is usually the case 


the preexistence of file QUAD is assumed. 
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e. Test Run 


To check whether the program is running correctly on 
any given system, Doleschall performs a limited calculation 
for n-D scattering using rank one singlet R-type potentials 
endeakuank onestensorpioreetonly. For comparison, he brings 
ehereourect resultsirombavcalculationegonehis home computer 
system. 

Tne CESeenungustearrtedeoureungrn ll by tnemcontrol tile 
JOBTEST which has the form of the JOB- files in the previous 
Section. OUEpuUttwEromttheselic cepied imtosthe permanent 
files QUADOUT, KEREOUT, FAMAOUT and PADEOUT and to *PRINT*. 

The run 1S restricted to one iteration of FAMA/PADE 
Cen J, 7=1/2 9ewithene concluding breakuptoupelasvic 
Scattering calculatrons.'’Becalisesitnisvattest rin, a larger 


ameunt of outputtrsecalled forvihantisenormaily required? 


f. Permanent Files 

There are a number of permanent files maintained in 
Gonnection wrth the running “Of thrs program. Tablesbaleiises 
those which have been stored on tape. data set refers to the 
FORTRAN data set referencersnmumber, if “any, with which the 
Pubes SeassSocrated dubingpexecuuion. 

In addition to the basic files of the program package, 
there are several files which relate specifically to the 
Dinynn)p reaction. Files preiaxed by OBS execuleva 
Simplified form of the program using mesh points from SQUAD 


and a rank one Yamaguchi potential. Otherwise they resemble 
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DOCODE 

DOLIBR 
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GAQUAD 

JOBBRUP 
JOBELAS 
JOBFAIT 
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JOBE ETS 
JOBKERE 
JOBNNFSI 
JOBNNOFS 
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description 


complete source listing 
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source listing of main programs 
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‘Symmetric Constant Relative Energy 
2special structure (see text) 
2simplified quadrature (using SQUAD) 
*modified’\n-p effective range 
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BhewuOBbriilesteitsection d, 

Files denoted by -RES- contain three-body T-matrices. 
BNese are associated with) file 2 of section ba They 
repmesent the result opecatculattonceup to the end of “the 
final FAMA-PADE iteration(19/2*" in all cases). For the more 
complicated potentials this involves a considerable amount 
of Gajeuieeie tenn making these files well worth retaining. 
Prefixes and suffixes indicate the potential used. 

NDRES contains three sets of T-matrices for each 
three-body angular momentum and parity: the first for the 
2°S 'R,4T4RA3,P,D interaceion,; eheesecond omitting D-waves, 
and the third omitting both P- and D-waves. This was 
constructed by defining the full potential in QUAD and KERE, 
anaeaunning PADE withStEsmtZ iow, Brive oe(three 
calculations). The potential weights (TOFI) were chosen so 
that P and D waves were omitted as required. The matrices 
are interleaved because FAMA and PADE iterations must be run 
in order of increasing angular momentum (JT). 

Physical resuleoeate produceds trom “RES- files “simply 
by sunning BRUP or ELAS @with the desired file on data Ser 2. 
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APPENDIX C PARAMETERS OF THE TWO-BODY INTERACTIONS 


Refer to table IV.2 for effective range parameters 


(S-wave only). 


The separable potential in the k,s(spin) partial wave 


has the form 
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The potential for coupling between k,s and m,t partial 


wave states has the general form: 
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